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In this paper, dedicated to the career of Iver Brevik, we review the derivation of the retarded van 
der Waals or Casimir- Polder interaction between polarizable molecules from the general multiple 
scattering formulation of Casimir interactions between bodies. We then apply this van der Waals 
QQ ' potential to examine the interaction between tenuous cylindrical bodies, including eccentric cylin- 

, ders, and revisit the vanishing self-energy of a tenuous dielectric cylinder. In each case, closed-form 

■ expressions are obtained. 
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' ^ '■ Since the earliest calculations of fluctuation forces between bodies [ij, that is, Casimir or quantum vacuum forces, 
multiple scattering methods have been employed. Rather belatedly, it has been realized that such methods could be 

1 used to obtain accurate numerical results in many cases 0, [1, H, 0] ■ These results allow us to transcend the limitations 

' of the proximity force theorem (PFT) 0,0] j so make better comparison with experiment, which typically involve 
, curved surfaces. (For a review of the experimental situation, see Ref. [1].) 

■ These improvements in technique were inspired in part by the development of the numerical Monte-Carlo worldline 
I— —I, method of Gies and Klingmiiller [l^, [HI but the difficulty with this latter method lies in the statistical 

limitations of Monte Carlo methods and in the complexity of incorporating electromagnetic boundary conditions. 

■ Optical path approximations have been studied extensively for many years, with considerable success [Tsl. ITsl. . 
However, there always remain uncertainties because of unknown errors in excluding diffractive effects. Direct numerical 
methods (TtI. ITsj. based on finite-difference engineering techniques, may have promise, but the requisite precision of 

. 3-dimensional calculations may prove challenging pJI]. 

' The multiple scattering formalism, which is in principle exact, dates back at least into the 1950s [20l[2l|. Particularly 
' I noteworthy is the seminal work of Balian and Duplantier (For more complete references see Ref. [13] ■) This 

I ' technique, which has been brought to a high state of perfection by Emig et al. [H], has concentrated on numerical 
] results for the Casimir forces between conducting and dielectric bodies such as spheres and cylinders. Recently, 
. we have noticed that the multiple-scattering method can yield exact, closed- form results for bodies that are weakly 
coupled to the quantum field [13] • This allows an exact assessment of the range of applicability of the PFT. 
The calculations there, however, as those in recent extensions of our methodology [2^, have been restricted to scalar 
' fields with ^-function potentials, so-called semitransparent bodies. (These are closely related to plasma shell models 
[li [Hi [13i [l^O The technique was recently extended to dielectric bodies [1^, characterized by a permittivity e. 
Strong coupling would mean a perfect metal, e — > oo, while weak coupling means that e is close to unity. 

In this paper we will give details of the formalism, and show how in weak coupling (dilute dielectrics) we recover the 
sum of Casimir-Polder or retarded van der Waals forces between atoms. Exact results have been found in the past in 
such summations, for example for the self-energy of a dilute dielectric sphere [13] or a dilute dielectric cylinder [3l| . 
Thus it is not surprising that exact results for the interaction of different dilute bodies can be obtained. It is only 
surprising that such results were not found much earlier. (We note that the additive approximation has been widely 
used in the past, for example, see Ref. [3^ . but here the method is exact. Also, there are many exact computations 
for non-retarded London forces betwen bodies, e.g., Ref. ^33i], but these results can only apply to very tiny objects on 
the nanometer scale.) In our previous letter [2^ we considered the force and torque on a slab of finite extent above an 
infinite plane, and the force between spheres and parallel cylinders outside each other. Here we will examine further 
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cylindrical geometries, such as concentric cylinders, and eccentric circular cylinders with parallel axes, in cases where 
the dielectric materials do not overlap. We will prove that the results can be obtained by analytic continuation of the 
energies found earlier for non-contained bodies. Finally, we will re-examine the self-energy of a dielectic cylinder [3ll |. 



II. GREEN'S DYADIC FORMALISM 

For electromagnetism, we can start from the formalism of Schwinger [s^l, which is based on the electric Green's 
dyadic F. This object can be identified as the one-loop vacuum expectation value of the correlation function of electric 
fields, 

r(r,i;r',0 = «(T{E(r,i)E(r',0}). (2.1) 

Alternatively, we regard the Green's dyadic as the propagator between a polarization source P and a phcnomcnological 
field E (where a:^ — (r, i)): 

E(a;) = J {dx')T{x, x) ■ P{x'). (2.2) 

We will only be contemplating static geometries, so it is convenient to consider a specific frequency as introduced 
through a Fourier transform, 

r(x,x') = J (*-*')r(r,r';c.), (2.3) 

in terms of which the Maxwell equations in a region where the permittivity e(u;) and the permeability fi{w) are 
constant in space read 

V X r = ioj^, V • * = 0, (2.4a) 

i V X * = -iujsT', V • r' = 0, (2.4b) 

where we have introduced F' = r + where the unit dyadic includes a spatial 6 function. The two Green's dyadics 
given here satisfy the following inhomogenous Helmholtz equations, 

(V^ + Lj^en)T' = X (V X 1), (2.5a) 

e 

(V^ +uj'^£H)^ = iujfiV X 1. (2.5b) 
In the following, it will prove more convenient to use, instead of Eq. (|2.5ap , 

^VxVx-£)F = l. (2.6) 

In the presence of a polarization source, the action is, in symbolic form, 

W^^Jp-T-P, (2.7) 

so if we consider the interaction between bodies characterized by particular values of e and /i, the change in the action 
due to moving those bodies is 

SW = ^ J P ■ ST -P = J E- SP-^ -E, (2.8) 
where the symbolic inverse dyadic, in the sense of Eq. (|2.6p . is 

F-i = X V X (2.9) 

that is, 

5F-F"^ = -F-^F~^ (2.10) 
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By comparing with the iterated source term in the vacuum-to- vacuum persistence amphtude exp iW, we see that an 
infinitesimal variation of the bodies resuhs in an effective source product, 



P(r)P(r') 

from which we deduce from Eq. ()2.7p that 



^isr-^, (2.11) 

3ff 



SW = -TrT • ST-^ = — Ti ST ■ = — STi InT, (2.12) 
2 2 2 ' ^ ' 

where the trace includes integration over space-time coordinates. We conclude, by ignoring an integration constant, 

W = -^Trlnr. (2.13) 

This is in precise analogy to the expression for scalar fields. Another derivation of this result is given in the Appendix. 
Incidentally, note that the first equality in Eq. (|2.12p implies for dielectric bodies {fi — 1) 

^^^"iJ^J idr)Se{r,to)Tkkir,r';Lo), (2.14) 
which is the starting point for the derivation of the Lifshitz formula [s^ in Ref. [33 |. 

III. REDERIVATION OF CASIMIR-POLDER FORMULA 

Henceforth, let us consider pure dielectrics, that is, set fi = 1. The free Green's dyadic, in the absence of dielectric 
bodies, satisfies the equation 



X V X -1 



■ — 



(3.1) 



so the equation satisfied by the full Green's dyadic is 

(r^-i - v)r = 1, (3.2) 

where V = e ~ 1 within the body. From this we deduce immediately that 

r = ii-ToV)-'ro. (3.3) 

From the trace-log formula p.l3p we see that the energy for a static situation {W = — J dtE) relative to the free-space 
background is 

E = ^Tr InTo 1 ■ T = -'-Tv ln(l - ToV). (3.4) 

The trace here is only over spatial coordinates. We will now consider the interaction between two bodies, with non- 
overlapping potentials, V — Vi + V2, where = Ea — 1 is confined to the interior of body a, a = 1, 2. Although it is 
straightforward to proceed to write the interaction between the bodies in terms of scattering operators, for our limited 
purposes here, we will simply treat the potentials as weak, and retain only the first, bilinear term in the interaction; 

E,2 ^ '-TrToViToV2. (3.5) 

Here, as may be verified by direct calculation [36], 

ro(r, r') = V X V X lGo(r - r') - 1 = ( VV - lC')Go(r - r'), (3.6) 
where the scalar Helmholtz Green's function which satisfies causal or Feynman boundary conditions is 

Go(r-r') = ^, i?=|r-r'|, (3.7) 
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the Fourier transform of the Euchdean Green's function, which obeys the differential equation 

(-V2 + C2)Go(r-r') = <5(r-r'), 

and C — —ibj. 

Thus the interaction between the two potentials is given by 



E\2 



1 [ dC 
2^ 



(dr)(dr') 



'^''47r|r-r'| 



V^i(r)V^2(r') 



The derivatives occurring here may be easily worked out: 

.-Kim 



+ \C\R + CR') + ^(3 + 3|C|i?| + CR') 



i?2 



ICI-R 



47ri?3 



and then contracting two such factors together gives 

\C\R 



(V,;Vj - C2<5y)^-^(V,V, - = (6 + I2t + lOt^ + + 2*4) 



(47ri?3)2 



where t = \C,\R- Inserting this into Eq, 

1 



we obtain for the integral over C, 



647r3i?7 Jo 



du e"" ( 6 + 6u + -u^ + -u^ + -u'* 
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(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



El2 — — 



r — r 



,/|7 



(3.13) 



which is the famous Casimir-Polder potential [37| . This formula is valid for bodies, which are presumed to be composed 
of material filling nonoverlapping volumes vi and V2^ respectively, characterized by dielectric constants ei and £2, both 
nearly unity. We emphasize that this formula is exact in the limit £1^2 — * 1, as discussed in Ref. (3^ . 



IV. ENERGY OF CYLINDER PARALLEL TO A PLANE 



In Ref. |29j we derived the energy of two uniform dilute cylinders, of radius a and h respectively, the parallel axes 
of which are separated by a distance R, R > a + b. In terms of the constant 



N = 



23 



6407r- 



■(£l-l)(£2-l). 



(4.1) 



the energy of interaction per unit length is 



327rA/' ^ " 2 ( "r^ 



■cyl-cyl 



1 f a^+b^ \ _ 1 f a^^b^ y 



i?6 



5/2 ■ 



(4.2) 



If we take R and b to infinity, such that Z = R — b is held fixed, we describe a cylinder of radius a parallel to a 
dielectric plane, where Z is the distance between the axis of the cylinder and the plane. That limit gives the simple 
result 



-yl-pl ^4 (i_„2/^2)5/2- 

This is to be compared to the corresponding result for a sphere of radius a a distance Z above a plane: 



sph— pi 



-N 



1 



z4 (i-ayz^y 



(4.3) 



(4.4) 



FIG. 1: Dielectric £2 hollowed out by a cylindrical cavity which contains an offset parallel dielectric cylinder ei. 



where v is the volume of the sphere. 

As an illustration of how the calculation is done, let us rederive this result directly from Eq. (|3.13p . We see 
immediately that the energy between an infinite halfspace (of permittivity ei) and a parallel slab (of permittivity £2) 
of area A and thickness dz separated by a distance z is 

dE ^dz , , ^, 

so the energy per length between the cylinder and the plane is 

^.r 2 , sine N-Ko^Z 

which is the result (14.31). 



V. ECCENTRIC CYLINDERS 

As a second illustration, consider two coaxial cylinders, of radii a and b, a < b. The inner cylinder is filled with 
material of permittivity ei, while the outer cylinder is the inner boundary of a region with permittivity £2 extending 
out to infinity. An easy calculation from the van der Waals interaction p.lSp 

64iv r p , , f^^ de 

Cco-cyl = — / dpp dp p 



3 7o ""A "Wo (p2 + p'2-2pp'cos0)3 

327rA^ , f°° , + Axy 

dx I dy- 



3 Jo Jb'^ {y-xf 



3(62 -a2)3- 



(5.1) 



This reduces to the dilute Lifshitz formula for the interaction between parallel plates if we take the limit ^ cx), 
a — > (X), with b — a — d held fixed: 

2NTrb EN , , 

€co-c.i--^, or (5.2) 

Note that the result (|5.ip may be obtained by analytically continuing the energy between two externally separated 
cylinders, given by Eq. (|4.2p . We simply take R to zero there, and choose the sign of the square root so that the 
energy is negative. That suggests that the same thing can be done to obtain the energy of interaction between two 
parallel cylinders, one inside the other, but whose axes are displaced by an offset R, with R + a < b, as shown in 

Fig.m 

16ttN (1 - a2/62)2 + (1 _^ ^2^52) j^2^j^2 „ 2Ryb'^ 

^ccc— cyl — o 7~i TTo" ■ (^"^J 

3 [(l-a2/52)2+ ^4/54 _ 2(1 + a2/62)^2/52]5/2 
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We can verify this is true by carrying out the integral from the Casimir-Polder formula (|3.13p 

on 7\r f^a f'OC /'27r /'27r 

€ccc-cyi = — 5— pdp p dp' / de dO' [p^ + - 2pp' cos{e -9') + R^- 2R{p cos 6 - p' cos 9')Y ■ 
•jtt Jo Jb Jo Jo 

.(5.4) 

The angular integrals can be done, but the remaining integrals are rather complicated. Therefore, let us simply 
expand immediately in the quantities x — pj p' and y — R/ p' ^ which are both less than one. Then we can carry out 
the four integrals term by term. In this way we find 

n.m— ^ ^ ^ ^ ^ ^ ^ ' 

which is exactly the series expansion of Eq. (|5.3p for small ajh and Rjh. 

By differentiating this energy with respect to the offset i?, we obtain the force of interaction between the inner 
cylinder and the outer one, 5^ = —d<t/dR. Evidently, that force is zero for coaxial cylinders, since that is a point 
of unstable equilibrium. For small _R, 5" grows linearly with R with a positive coefficient. The inner cylinder is 
attracted to the closest point of the opposite cylinder. Similar considerations for conducting cylinders were given 
in Refs. [ssl [39|. with the idea that the cylindrical geometry might prove to be a useful proving ground for Casimir 
experiments. 



VI. SELF-ENERGY OF DILUTE CYLINDER 

This is a rederivation of the result found by dimensional continuation in Ref . |3l| . The summation of the Casimir- 
Polder forces between the molecules in a single cylinder of radius a is given by (in A'^ £1 = £2) 

327V r , r , / / de 

^'^'^-^ Jo '''Jo''' Jo {p^ + p-~2pp'cosef 

32 / du + (6.1) 



^0 



This is, of course, terribly divergent. We can regulate it by analytic continuation: replace the highest power of 
(p2 — p'^)^^ — {xu)~^, 5, by (3, and regard /3 as less than 1. Then, 



p _ 16iV 







dxx'-^ / duiu^-f" -6u'-^ + 6u-^) 



_ _16N 2,4-g (5-/3) 

3 (i_/3)(2-/3)(3-/?)- ^^-"^ 

Now if we analytically continue to /3 = 5 we get an vanishing self-energy to order {e — 1)^. This result was first 
discovered by Romeo (private communication) , verified in Ref. [3l| , and confirmed later by full Casimir calculations 
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APPENDIX A: DERIVATION OF GREEN'S DYADIC FORMALISM FROM CANONICAL THEORY 



Here we begin by sketching the development of the Green's dyadic equation from canonical quantum electrodynam- 
ics. For simplicity of the discussion, we will consider a medium without dispersion, so that s and p are constant. First 



we must state the canonical equal-time commutation relations. We will require (only transverse fields are relevant) 

[E,{v,t),E,{v',t)]^0. (Al) 

In a medium, it is the electric displacement field which is canonically conjugate to the vector potential, so we have 
the equal-time commutation relation (Coulomb gauge) 



[A(r, t), aoA(r', t)] = \U- ^ ) '5(r - r'). 



Now in view of Eq. (|2.ip , and the Maxwell equations 



(A2) 



V X E 



dt 



B, 



V X is = ^eE, 
^ ot 



we deduce from Eq. (|2.1 



(A3a) 
(A3b) 



V X V X r' + en—V = le^iSit - t')([E(r, t), E(r', t)]). 



But according to Maxwell's equations and Eq. (|A2 



(A4) 



[E(r, t), E(r', t)] —V X V X [A(r, <), -aoA(r', t)] = — ^ V X V X l^(r - r') 



(A5) 



If we now insert this into Eq. (jA4p we obtain for the Fourier transform of the Green's dyadic 



(V X V X -e/iw2)r' = iv X V X 1, 



(A6) 



which is indeed the equation satisfied by the solenoidal Green's dyadic, Eq. (I2.5ap . 

It is equally easy to derive the trace-log formula. We have the variational statement, for infinitesimal changes in 
the permittivity and the permeability [43 |. 



SE^-^ J {dr)(SeE'^ + SfiH^). 



(A7) 



Given Eqs. and (|Ma]) we can write this in terms of the coincident-point limit of Green's dyadic, 



SE=- I (dr) 



Se- S { - ] -^V X Vx 



r(r,r') 



(A8) 



according to Eq. (|2.9p . which involves an integration by parts, and coincides with the first equality in Eq. (|2.12|1 . 
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